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Low field vortex dynamics over seven time decades in a Bi 2 Sr 2 CaCu2 08+5 single 

crystal for temperatures 13 K < T < 83 K 

M. Niderost, A. SuteiE, P. Visani, and A.C. Mota 
Laboratorium fur Festkorperphysik, Eidgendssische Technische Hochschule Zurich, 

8093 Zurich, Switzerland 

G. Blatter 

Theoretische Physik, Eidgendssische Technische Hochschule Zurich, 
8093 Zurich, Switzerland 
(February 1, 2008) 

Using a custom made dc-SQUID magnetometer, we have measured the time relaxation of the 
remanent magnetization M rem of a E^S^CaC^Os+s single crystal from the fully critical state 
for temperatures 13 K < T < 83 K. The measurements cover a time window of seven decades 
10~ 2 s < t < 10 s, so that the current density j can be studied from values very close to j c 
down to values considerably smaller than j c . From the data we have obtained: (i) the flux creep 
activation barriers U as a function of current density j, (ii) the current- voltage characteristics 
E(j) in a typical range of 10 -7 V/cm to 10~ 15 V/cm, and (hi) the critical current density j c (0) at 
T = 0. Three different regimes of vortex dynamics are observed: For temperatures T < 20 K the 
activation barrier U(j) is logarithmic, no unique functional dependence U(j) could be found for 
the intermediate temperature interval 20 K < T < 40 K, and finally for T > 40 K the activation 
barrier U(j) follows a power-law behavior with an exponent jj, ~ 0.6. From the analysis of the data 
within the weak collective pinning theory for strongly layered superconductors, it is argued that for 
temperatures T < 20 K pancake- vortices are pinned individually, while for temperatures T > 40 K 
pinning involves large collectively pinned vortex bundles. A description of the vortex dynamics in 
the intermediate temperature interval 20 K < T < 40 K is given on the basis of a qualitative low 
field phase diagram of the vortex state in Bi2Sr2CaCu20g+i. Within this description a second peak 
in the magnetization loop should occur for temperatures between 20 K and 40 K, as it has been 
observed in several magnetization measurements in the literature. 

PACS numbers: 74.25.Dw, 74.60.Ge, 74.60.Jg, 74.72.Hs 



I. INTRODUCTION 

High temperature superconductors (HTSC) are char- 
acterized by large values of the Ginzburg-Landau param- 
eter k = A/£, so that most of the H — T phase diagram is 
dominated by the presence of vortices. Furthermore, the 
high anisotropy of P^S^CaC^Os+a (Bi2212) has strong 
implications for the behavior of the flux lattice in the 
mixed state. When a magnetic field H is applied perpen- 
dicularly to the afr-planes, the vortices.-can be described 
as two dimensional "pancake-vortices"B lying in the su- 
perconducting Cu02 layers. These pancake-vortices in- 
teract both through the interlayer Josephson coupling as 
well as through electromagnetic coupling. Such a layered 
vortex structure is very sensitive to thermal and quantum 
fluctuations, especially considering the small coherence 
length £ in the direction parallel to the Cu02 planes. As 
a consequence, pinning is relatively weak as compared to 
classical type II superconductors and strong relaxations 
of the magnetization M are observedtTEl which deviate 
from a pure logarithmic time dependence. 

Since the discovery of the HTSCo, the theoretical and 
experimental work concerning vortices and their dynam- 
ics have been strongly intensified. Those investigations 
were mainly focused on a regime where the current den- 



sity j is relatively small as compared-top-the critical cur- 
rent density j c . Only little is knownt2rt3 at present re- 
garding the vortex dynamics in a regime where the cur- 
rent density j is close to j c . 

In the here presented work, we investigate experimen- 
tally the low field vortex dynamics in a Bi2212 single crys- 
tal for magnetic fields H perpendicular to the Cu02 lay- 
ers. For this purpose we have designed and constructed a 
dc-SQUID magnetometer with high sensitivity and long 
time thermal stability. The measurements of the relax- 
ation of the remanent magnetization M rcm are taken in 
the temperature interval 13 K < T < 83 K and cover a 
time window of seven decades. The wide current range 
of the experimental data allows a detailed analysis of 
the vor tex d ynamics within the theoretical vortex creep 
models.liiali-3 Applying the method of Maley et al&3 to the 
relaxation data, a characteristic functional dependence 
between the activation barrier U and the current den- 
sity j is obtained for the temperature regimes T < 20 K 
and T > 40 K, whereas for temperatures between 20 K 
and 40 K the [/(j)-relation is found to depend strongly 
on temperature.. On the basis of a qualitative low field 
phase diagramtEl of Bi2212, an interpretation of the be- 
havior of the vortex dynamics in tie temperature interval 
around the crossover temperaturell3i23 T ~ 25 K is given. 
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II. EXPERIMENTAL 

The measured single crystal is 0.9 x 1.3 x 0.05 mm 3 in 
size and the critical temperature T c is 95 K, as deter- 
mined by ac-susceptibility measurements. The growth 
procedure as well asr-fche transport properties have been 
described elsewhere.EH The experiments are performed 
in a custom made dc-SQUID magnetometer, where the 
sample remains stationary in the pick-up coil during the 
measurements. The externally applied magnetic field H 
is supplied by a superconducting coil, working in a non- 
persistent mode. In order to prevent eddy currents, the 
experimental cell is entirely built out of epoxy resin Sty- 
cast 1266E3 

During the measuring procedure, the sample is first 
zero field cooled in the residual field of the cryostat 
(H TCS ~ lOmOe) from well above T c and then stabi- 
lized at a fixed temperature T. Next, a magnetic field 
H applied perpendicularly to the a6-planes is gradually 
increased from zero to iJ max before being removed lin- 
early at a rate of 9T/s. This fast rate is achieved by 
shorting the superconducting coil [L ~ 7H) over an ex- 
tremely non-linear resistor. Measurements of the current 
in the coil show that no discontinuities occur during the 
removal of the field. The data are taken as soon as the 
decreasing magnetic field fulfills the condition H < 1 Oc. 
As a time origin for the measured data we choose the 
time at which H starts being removed. After measuring 
the relaxation of the remanent magnetization M rcm for 
about seven time decades, the sample is heated above T c 
in order to record its residual magnetization. The maxi- 
mum values of the initially applied magnetic field 7J max 
are shown in Table [| for all the measuring temperatures 
T. The values of if max are chosen so that they are bigger 
than twice the field needed to achieve full flux penetra- 
tion into the sample. Table [| further contains the values 
of the residual field H rcs along the axis of the supercon- 
ducting coil after cycling the magnetic field from zero to 
iJ max and back to zero again. 

Due to the high field removal rate H = 9T/s, it was 
necessary to perform some controls concerning the initial 
field profile in the sample as well as self-heating effects. 
Several field removal rates H have been tested. We found 
that for the field removal rates 10~ 2 T/s < H < 9T/s, 
the measured remanent magnetizations M rcm do not 
show any remarkable difference. Moreover, no significant 
change in the dynamics of the relaxation of M rem could 
be detected by increasing the initial field values iJ max by 
a factor of 2-3. From our estimations we concluded that 
the self-heating of the sample due to induction as well 
as to flux flow can be neglected for all temperatures and 
fields of our measurements. 

With the described experimental procedure we mea- 
sured the relaxation of the remanent magnetization M rem 
for the Bi2212 crystal from the fully critical state in a 
temperature range of 13 K < T < 83 K. Since there is 
only a small uncertainty of the time origin of the creep 



process (< 18 ■ 10~ 3 s), the initial behavior of the relax- 
ation data as a function of time is very well defined so 
that we were able to test the existing vortex creep models 
over a wide current density region starting from values 
near j c . 



III. FLUX DYNAMICS MODELS 

The main effect of pinning is to allow a flux density 
gradient to be sustained within a type II superconduc- 
tor. This is intrinsically related to the flow of a macro- 
scopic diamagnetic screening current density j that can 
be expressed, in the continuous approximation, through 
Maxwell's equation VAB = 4n/c j. The configura- 
tion with a finite flux density gradient is metastable and 
hence is bound to decay. The dynamics arises from the. 
vortex creep motion as.-a. .result of thermal activationEj 
and quantum tunnelingt-xEj (T < 5K). For a geometry 
where B \\ z and j \\ y, the Maxwell equations together 
with the condition of flux. .conservation lead to the non- 
linear diffusion equationcSEll 



Of 



47T dx 2 



(vB). 



(1) 



Andersoncj postulated that the velocity of the vortices, 
as a consequence of thermal activation over the pinning 
barrier U(j), be given by 



(2) 



where v a (j) is the mean vortex velocity and can be ex- 
pressed as v {j) = l{j)/T , where is the mean hop- 
ping length and r the is inverse attempt frequency. For 
a situation where v {j) is independent of j, the diffusion 
equation ([!]) can be transformed into 



2i 
at 



Jc 

— exp 



urn 

k B Tj 



(3) 



As discussed by Geshkenbein and LarkinH, equation (j^) 
can be solved within logarithmic accuracy, yielding 

TABLE I. Maximum applied magnetic field # max for dif- 
ferent measuring temperatures T and values of the residual 
field H ICS , which is due to the flux remaining trapped in the 
superconducting coil after the removal of the external field H . 
The residual field of the cryostat is about 10 mOe in opposite 
direction to the applied magnetic field. 



T(K) 


13 


15 - 27 


30 - 40 


50 


> 60 


# max (Oe) 


1500 


1000 


500 


300 


< 100 


H ies (mOe) 


710 ± 50 


480 ± 50 


60 ± 10 


20 ± 10 


-10 ± 10 
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U(j(t))~k B Tln(l + 



(4) 



where t = ksTT / 'j c \djU\ is a time scaling factor. Once 
the functional dependence between the pinning barrier U 
and the current density j is known, the time dependence 
of j is simply determined by the inversion of (^) . 

On approaching the critical current density j c , the ef- 
fective pinning barrier vanishes and one can write 



U(j -> j c ) ~ U c ^1 

Comparing equations (Q) and (Jj 
pendence of j is obtained: 



(5) 



the following time de- 



3(t) Jc 



In 1 



l/a 



j -> Jc (6) 



for 



which maps to the original formulation of Anderso 
a = 1. 

In the above derivation it is assumed that the current 
densities j be close to j c . This is a good assumption for 
conventional type II superconductors. Further theoreti- 
cal considerations are necessary to describe the strongly 
decaying current densities in HTSC, for which values of 
j much smaller than j c are reached already at laboratory 
times. For the HTSC in the lirnit of small currents, the 
weak collective pinning theoryo (WCPT) as well as the 
vortex glass theoryO, predict an activation barrier that 
diverges algebraically for vanishing currents: 



' (y 



(7) 



Inserting the relation (Q) into equation (Q) the following 
non purely logarithmic time dependence of the current 
density j is obtained: 



]{t)-Jc 



1 -1/M 



j < jo- 



(8) 



In order to find a more general formula, (g) and (g) 
(we assume a = 1) can be interpolated with the follow- 
ing expression 



j{t)-jc 



1 



(9) 



and the corresponding activation barrier is (see equation 



)) 



J - -1 



(10) 



Within the single vortex pinning regime the exponent 
1/fX w 7 is large, such that for fj, ksT/Uc In (1 + t/t Q ) <C 
1 expression (O) can be approximated by 



j(t) ~ j c 1 + 



f 

with a logarithmic potential 



-k B T/Ua 



U(j)~U c ln(j c /j). 



(11) 



(12) 



Notice that within the WCPT, the divergence in the 
potentials (|?|), (pTj|), and ( ^2[ ) at low current densities j 
is related to the observation that the activated motion of 
vortices involves hops of larger vortex segments /bundles 
over longer distances. The elastic energy cost will there- 
fore grow with decreasing j. This is no longer the case 
for the point-like pancake- vortices for which no extra de- 
formation energy is needed in order to overcome the pin- 
ning barrier for decreasing current densities. For strongly 
layered superconductors within the single pancake creep 
regime the activation barrier U(J) is therefore expected 
to saturate. .However, using the concept aLvariable-range 
hoppingE3E!l (VRH) it has been arguedES that, for de- 
creasing current densities j, pancake- vortices still couple 
into a 2D elastic manifold. As a matter of fact, due to the 
randomness in the energies of the mctastable state, pan- 
cakes will hop over larger distances as the current density 
j decreases. Such a large hopping distance u leads to a 
large shear interaction energy cqqoIu 1 between pancakes 
(c66 is the shear modulus and d the interlayer distance). 
As a consequence, for a large enough hopping distance it, 
the pancake- vortex will start to couple to its neighbours. 
Thus, the vortex system is expected to first go through a 
VRH ijegime, which is followed by a 2D collective creep 
regimeQ at still lower current densities. 

In the above treatment of the flux dynamics models 
we have considered current densities j flowing inside a 
superconductor, whereas from the experiment we obtain 
spatially averaged values of the magnetization M. In 
the case of an infinite slab parallel to the applied mag- 
netic field H, the dependence between M and j has been, 
described by BeanEj Recently Gurevich and BrandtEll 
obtained an asymptotic solution for the non-linear dif- 
fusion equation |lj) describing flux creep in strips and 
disks starting from a barrier as given in formula ([l0]). It 
turns out that, despite the particular field distribution 
for these sample geometries, the current density j can 
still be considered as constant throughout the sample at 
a given time t. It follows that the magnetization M, 
which is given by 



M(t) = ~y c / Jlr.D.IV. 



(13) 



for a disk-like geometry and for a constant current den- 
sity j(r, t) = j(t) e$, can be expressed as 



|M(i)|=j(i)-i-i J \rAet\dV, 



(14) 



where the integration over the geometrical factor leads 
to 



3 



M(*)^-j(i), 



(15) 



with R being the sample radius. For the case of disks 
(strips) the well known Bean model relationship for an 
infinite cylinder (infinite slab) in the fully critical state 
is therefore still a valid approximation. 



IV. VORTICES IN STRONGLY LAYERED 
SUPERCONDUCTORS 

For the considerations given in this Section concern- 
ing the vortex lattice in coupled superconducting lay- 
ers we will closely follow the approach of Refs. |l3| and 
|34]. Within weak collective pinning theory the size of the 
correlated regions (Larkin domains) is determined by the 
balance between deformation energy and pinning energy. 
In terms of length scales, the volume forming the Larkin 
domain is given by the pinning correlation lengths R c and 
L c in the direction perpendicular and parallel to the mag- 
netic field, respectively. Through the study of the relative 
magnitude of the deformation and the pinning energy of 
a vortex lattice in coupled superconducting layers, it is 
possible to determine the size of the correlated regions as 
a function of temperature and field. 

For a magnetic field H perpendicular to the su- 
perconducting layers, a vortex lattice has three rel- 
evant energy scales, namely the tilt energy [/tilt ~ 
C44(R C ) ■ (r p /L c ) 2 R 2 L c , the shear energy [/ shear ~ 
C66 • (r p /R c ) 2 R 2 L c , and the pinning energy U pin w 
{^^R^Lc/rpCL 2 ) 1 / 2 , with c 44 being the dispersive tilt 
modulus, C66 the shear modulus, r p (T) the range of the 
pinning force, a the intervortex spacing, and 7 the disor- 
der strength (where a short-scale correlated disorder po- 
tential has been assumed (U p i a (r), U p i n (r')) = 76(1 — r')). 
Depending on the relative magnitude of these energies, 
one can distinguish four possible pinning regimes: (1) in- 
dependently pinned vortex pancakes (OD pinning regime: 



Upm > Unit, [/shear), (2) independently pinned vortex 
lines (ID pinning regime: U t ilt > U pin > [/shear), (3) a 2D 
collectively pinned state in which the 2D vortex lattices 
in the layers are pinned independently from each other 
(C/shear > U p i n > Unit), and (4) a 3D collectively pinned 
state (U t at, U shcar > U pm 

According to 
T < 



Refs. 

j2 , E n 



13 



and | 
and 



lq , for temperatures 
T m {U' 2 c E pc y^ and fields B < B 02 « 
10 $ /(27t£ 2 ) (j c (0)/ja), where T is a few tens of Kelvins 
and Bq2 is a few Teslas, the dominant energy scale 
for the strongly anisotropic Bi2212 is the pinning en- 
ergy Upin (where U pc ~ e d(j c /j ), E pc » e d(£/X) 2 , 
£o = ( < &o/47rA) 2 , <i> = hc/2e, d is the interlayer distance, 
and j the depairing current density). The B-T phase 
diagram for this region is therefore characterized by OD 
pinning. On the other hand, for temperatures T > Tq the 
collective pinning length L c and the collective pinning ra- 
dius R c both grow very fast due to thermal depinning. 
This implies that for temperatures T > 20 K the size of 



the Larkin domains becomes large giving rise to a 3D 
pinning regimeEl At high fields B > B 2 3,r^ crossover 
to a 2D collective pinning region is predictedl!3 when the 
shear energy outweighs the tilt energy. 

Finally, since the relaxation measurements of the re- 
manent magnetization M rom presented in this work are 
performed in the "field off" state, we need to discuss the 
very low field regime. At fields B > <I> /A 2 , the shear 
modulus C66 has a linear dependence in B, wherq a^ , a, t 
low fields (B < $ D /A 2 ), cqq decreases exponential! 

2\ 1/4 



C66 ~ < 



|(^) e-V^, B<$ o/ A 2 , 



e B 
I 4$ ' 



B > $ /A 2 , 
(16) 



where A is the penetration depth. As a consequence, also 
the shear energy C/ S hear decreases exponentially for fields 
B < $ Q /A 2 . This means that for temperatures T > Tq 
and small enough magnetic fields (B < Bi 3 ) a ID pin- 
ning regime occurs. Fig.[l] shows a qualitative map of the 
low field pinning regimes of Bi2212 resulting from these 
considerations. 




FIG. 1. Qualitative low field phase diagram of the vortex 
state in Bi2212 for magnetic fields H perpendicular to the 
superconducting layers. The differently shaded areas in the 
Figure represent the following pinning regimes: OD, individu- 
ally pinned pancake- vortices; ID, individually pinned vortex 
lines; 2D, collectively pinned state in which 2D lattices of 
pancake- vortices in the layers are pinned independently from 
each other; 3D, collectively pinned vortex bundles. B13 (B23) 
represent the fields at which the ID (2D) regime crosses over 
to the 3D regime. To is a crossover temperature terminating 
single pancake pinning. A sketch of the irreversibility line IL 
and of the upper critical field H C 2 is also given. 
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FIG. 2. Normalized remanent magnetization vs. time, 
measured after cycling the sample in an external magnetic 
field H . In parenthesis the values of the maximal cycling 
fields # max are given: 

O (tf max = 1 kOe) at 15 K, □ (# max = 1 kOe) at 25 K, and 
O (# max = 40 Oe) at 69 K. 

The time origin is given by the instant, when the exter- 
nally applied magnetic field H starts being decreased and 
t s — 18 • 10 -3 s is the time, when the first point of the relax- 
ation of M rcm is taken. 



analysis of the time evolution of the current density j is 
given. 

A. Low Temperature Regime (T < 20 K) 

As shown by Maley et al. it is possible to determine the 
activation barrier for vortex motion U(J) directly from 
the relaxation data j(t). Starting from equation (0), one 
obtains 



U(j) 



-IcbT In 



dt 



ksT In 



(17) 



where the term ksT ln|sj c /r | is independent of j, and 
s = 1 cm 2 s/A. Plotting the expression -fc^T In \s dj/dt\ 
as a function of current density at different temperatures 
T, a set of curves is found which are vertically shifted 
with respect to each other. For a temperature interval 
where the functional dependence between the activation 
barrier U and the current density j is essentially temper- 
ature independent, this shift is given by the term aAT, 
where a ~ ln|sj c /r | is a constant, and AT = Ta — T\ 
is the temperature difference between two considered 
curves. Combining the data measured at different tem- 
peratures T, the activation barrier U(j) is obtained over 
a wide current density range. 



V. EXPERIMENTAL RESULTS AND ANALYSIS 

In Fig. U we illustrate the time dependence of the re- 
manent magnetization A/ rcirr -.with a typical set of data. 
A non-logarithmic behaviorE3 is observed at all temper- 
atures. Notice that at 25 K, where a sharp drop in the 
relaxation rate—Sp- ■= — d\nM Tcm /dlnt has been previ- 
ously reportedtZrt3, the remanent magnetization M rem 
decays extremely fast in the first few seconds after the 
removal of the external field H. This is also seen in 
Fig. [|, where the current density j (as obtained from for- 
mula (|l5|)) is plotted as a function of temperature for the 
times t s ~ 18 • 10~ 3 s, h = 1 s and t 2 = 10 4 s. The data 
taken at the starting time t s ~ 18 • 10~ 3 s (empty cir- 
cles) suggest the presence of only two regimes of vortex 
dynamics, separated by a crossover at T ~ 30 K. Both 
regimes are accurately described by an exponential tem- 
perature dependence but with different slopes dhij/dT. 
However, at longer times t > 1 s (filled circles and empty 
diamonds) the existence of a third regime for temper- 
atures between 20 K and 40 K becomes evident. This 
third regime is characterized by very particular vortex 
dynamics and will be referred to as the "intermediate 
regime" . We will discuss these temperature regimes sep- 
arately and distinguish them as follows: a low temper- 
ature regime for T < 20 K, an intermediate regime for 
20 K < T < 40 K and a high temperature regime for 
T > 40 K. For each regime we determine the activation 
barrier U(j) by means of the method of Maley et alx3 
Once the functional dependence of U(j) is obtained, an 
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FIG. 3. Current density j as a function of temperature for 
different times t: 

O starting time t s ~ 18 ■ 10" 3 s, • ti = Is, and O t 2 = 10 4 s. 
The lines serve as guides to the eyes. 

For temperatures T < 20 K, the data obtained from 
the expression — k B T In \sdj/dt\ at different tempera- 
tures T can be accurately mapped onto a common curve 
using a single constant a. The obtained potential U (j) is 
shown in Fig.^. It is interesting to observe in Fig. f|(a), 
that the data measured at a fixed temperature T (marked 
by horizontal segments) do overlap over wide regions of 
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current. As seen in Fig.|(b), the potential U(j) is pro- 
portional to the logarithm of the current density j over 
a wide current region. This is in good agreement with, 
previous relaxation measurements by van der Beek et alH 
and by Emmen et alu, who found a logarithmic depen- 
dence of U(j) for temperatures 4K < T < 17K. The 
deviation from the logarithmic behavior at temperatures 
T > 19 K is attributed to the influence of the approaching 
intermediate regime. 
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FIG. 4. (a) Flux creep activation barrier for temperatures 
13 K < T < 23 K as determined from the magnetic relaxation 
data by the method of Maley et al. (the constant used for 
matching the curves is a = 26 ± 1). The horizontal segments 
represent the current windows as obtained from the data at 
a fixed temperature T. 

(b) The same data in a semi-logarithmic graph. The line is a 
fit for temperatures up to T — 19 K (indicated by the arrow) 
with a logarithmic potential of the type U(j) — U c la(j c / j). 
From the fit one finds the T = critical current density 
j c (0) ~ (1.0 ± 0.3) ■ 10 6 A/cm 2 . 

For temperatures T < 19 K, a fit of the measured 
potential U(J) with the logarithmic activation barrier 
( |j~2| ) leads to the following parameters: U c ~ 140 K and 
the extrapolated critical current density j e (T = 0) ~ 
1 • 10 6 A/cm 2 . The value of jJT = 0) is very close to the 
values found in the literature&Q (taking into account the 
considered proportionality factors between M and j). 



As discussed in Section HI, the activation barrier U(j) 
is expected to be logarithmic within the single vortex 
pinning regime. Since the measured potential U(j) is 



indeed logarithmic, this would suggest that for temper- 
atures T < 20 K vortex strings are pinned individu- 
ally. However, a simple estimate of the collective pin- 
ning length along the c-axis L° c ~ £^(jo/ jc) 1 ^ 2 , where 
e is the anisotropy factor and j is the depairing cur- 
rent density, shows that, for the parameters of Table || 
and H M ~ 2 A < d = 15 A. This means that, for 



temperatures T < 20 K and low enough magnetic fields, 
pancake- vortices placed on different superconducting lay- 
ers are pinned independently indicating the presence of 
a single pancake pinning regime. A more detailed dis- 
cussion of the low temperature activation barrier will be 



given in Section VI 



We can crosscheck the result for the barrier as obtained 
via the Maley analysis making use of equations ( pd| ) and 
(|l5| ) . A typical fit to the data measured at temperatures 
T < 19 K is shown in Fig.|5|, confirming the logarithmic 
dependence U(j) ~ U c ln(j c /j). The resulting fitting pa- 
rameters are the following: U c — 140 K, t Q ~ 3 • 10 -2 s 
and values of j c (T) about 5% above the values shown in 
Fig. | for t s ~ 18- 10~ 3 s. 

Finally, we point out that for temperatures T < 20 K 
the values of the pinning potential U(j) and of the ex- 
trapolated critical current density j c (T — 0) are both in 
good agreement with the results in the literature, usu- 
ally obtained in the "field on" mode at much slower field 
ramping rates H. 
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FIG. 5. Normalized remanent magnetization vs. time for 
temperatures T = 17 K and T = 69 K. The lines are fits ac- 
cording to formula ([H]) for the 17 K data and formula Q for 
the 69 K data. 



B. Intermediate Regime (20 K < T < 40 K) 

In order to find the activation barrier U(j) for temper- 
atures 20 K < T < 40 K, the relaxation data are again 
evaluated with the method of Maley et al. The results 
obtained with help of equation ( |l7| ) for different temper- 
atures T are shown in Fig. |[ We observe that the curves 



G 



are strongly tilted with respect to each other and it is not 
possible to obtain a unique smooth curve by simply shift- 
ing the data obtained at different temperatures T along 
the vertical axis. Thus, within the temperature range 
20 K < T < 40 K, we cannot find a unique temperature 
independent functional relation between U and j follow- 
ing the above approach. A qualitative interpretation of 
the vortex dynamics in this temperature regime will be 



given in Section VI 
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FIG. 6. Flux creep activation barrier vs. current density j 
for temperatures 23 K < T < 35 K. The vertical axis is only 
defined up to a constant value. For this temperature regime 
the curves are strongly tilted with respect to each other and 
cannot be "glued" onto a common curve anymore. 



C. High Temperature Regime (40 K < T < 83 K) 

For temperatures T > 40 K, the activation barrier U (j) 
is found with the same method that has been applied for 
the low temperature regime using a single constant a. 
The resulting barrier U(j) is shown in Fig.^. From the 
double logarithmic plot of Fig. 0(b), we observe that the 
activation barrier U(j) follows a power-law behavior over 
a wide current range. Fitting this potential with formula 
@ for temperatures 62 K < T < 83 K, we find the val- 
ues U c ~ 1000 K and fi ~ 0.6. 

Accor ding to Ref. [l3|, a power-law potential with the 
form of ([n]) leads to a time dependence of the current 
density j as given by the interpolation formula (^J). As 
one can see from the solid line of the fit to the T = 69 K 
data in Fig.|EJ the time dependence of the current density 
j is very well described by the interpolation formula (|^). 
The fitting parameters confirm the results previously ob- 
tained for the barrier and can be summarized as follows 
(see also Table 0): U c ~ 1000 K, fi ~ 0.6, t ~ 3 • 10~ 3 s, 
and values of j c (T) about 5% above those obtained from 
Fig. H at t s ~ 18 • 10~ 3 s. According to weak collective 
pinning theory, an exponent /i ~ 0.6 indicates a regime 
of large 3D bundle pinning. 

The high temperature data have been analysed consid- 
ering a constant current density j inside the sample 
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FIG. 7. (a) Flux creep activation bar- 
rier for 40 K < T < 83 K as determined by the method of 
Maley el al. (with a = 62 ± 1) from the magnetic relaxation 
data. The horizontal segments represent the current windows 
as obtained from the data at a fixed temperature T. 
(b) The same data in a double-logarithmic graph. The line is 
a fit for temperatures T > 62 K, with a power-law potential 
as given in formula (p"o[). 



as assumed in the Bean model (see formula (p.5|)). We 
argue that, for the present measurements of the relax- 
ation of the remanent magnetization M rem , the contri- 
butions of pinpmg-.due to potential barric^ .arising from 
surface effectsE3~Ea and sample geometryEl'liil have only 
a secondary effect as compared to the contributions of 
bulk pinning. As a matter of fact, if surface barriers 
were the only mechanism responsible for the irreversible 
behavior, the magnetization curves would be character- 
ized by zero-jjiagnetization on the descending branch of 
the loopO'LJ Fig.|| shows three magnetization cycles of 
the Bi2212 crystal measured at different temperatures T. 
From the shape of the curves we can safely say that, for 
our sample, pinning due to surface barriers does not play 

TABLE II. Experimental fitting parameters obtained from 
the relaxation data in the low temperature regime (T < 20 K) 
and in the high temperature regime (T > 40 K). 





Uc (K) 


t (s) j c (T = 0) (A/cm 2 ) 




T < 20 (K) 


140 


3 • W~" 1 ■ 10 b 





T > 40 (K) 


1000 


3- 10" 3 


0.6 
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a dominant role. In thin superconducting strips of rect- 
angular cross sectiap^ Meissner currents flow throughout 
the whole samplec3c3 and not only in a surface layer of 
width A. Lorentz forces arising from the Meissner cur- 
rents will therefore act on the vortices and give rise to a 
barrier of purely geometric originO This kind of geomet- 
rical barrier will not influence measurements performed 
in the "field off" state, since Meissner currents do not 
flow in this state (the influence of the residual fields H res 
is of minor importance, see Table |). 



200 
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FIG. 8. Magnetization curves at different temperatures T 
for the Bi2212 crystal. The descending branches of the loops 
do not have zero magnetization so that we conclude that 
Bean-Livingston surface barriers are only of secondary im- 
portance. 



D. Current- Voltage Characteristics 

Inductive measurements of the relaxation of the, Rema- 
nent magnetization M rcm are a powerful toolH'Eacia for 
the evaluation of E{j) characteristics of a superconduc- 
tor down to very low values of the electric field E. The 
functional dependence between the electric field E and 
the current density j can be found as follows. In the 
Coulomb gauge V • A = 0, the vector-potential A can 
be expressed as 



A(r) 



'3,,.' 



d A r 



where j is the current density. For a disk-like geometry 
and for a constant current density j(r') = j e^, (j = 
const.), one obtains 



A(r) 



1 

J ■ - 

c j v 



(19) 



which, together with the Faraday induction law dtA 
—cE, leads to 
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FIG. 9. Current-voltage characteristics as extracted from 
the relaxation data for different temperatures T at zero ex- 
ternal field H . 



E = 



c z at 



7 



(20) 



where E is the mean electrical field in the sample, R the 
sample radius and h its thickness. The constant factor 
7 arises from the assumption of a disk-like geometry for 
the sample. 

The results obtained from the magnetic relaxation data 
and equation (|2(]) are plotted in Fig. [)[ From the graph 
one can clearly distinguish the three regimes of vortex 
dynamics, which were previously discussed in this Sec- 
tion. For current densities j < j c , the electric field E due 
to the dissipative process of thermally activated vortex 
drift can be written asa 



E(j) = E c exp [-U(j)/k B T] 



(21) 



From Fig.|| it is seen that for temperatures T < 19 K 
the E(j) curves follow a power-law behavior. As a mat- 
ter of fact, inserting the logarithmic potential ( |T^ ) into 
formula (|2l]) one obtains 



E 
E~, 



U c /k B T 



(22) 



From a fit to the curves in this temperature regime we 
obtain again U c ~ 140 K, in agreement with the previous 
results. 

TABLE III. Values of j„, j o /j c (0), L c c (0), E pc , and U pc as 
(18) obtained from following formulas: j a — c<I> /(12v / 37r 2 A 2 ^), 
L c c =a eajo/jc) 1/2 , E pc » e d(C/A) 2 , and U pc ~ e d(j c /j ). 
The parameters used for the theoretical estimates are given 
for the configuration where H is perpendicular to the super- 
conducting layers. 



jo w 10 8 (A/cm 2 ) U pc ~ 20 (K) L c c (0) ~ 2 (A) 

j o /j c (0)^100 E pc « 3 • 10- 2 (K) 
Xl — 1800 (A) 
^bcs =i 30 (A) 



A(0) = \ L /V2^ 1300 (A) d ~ 15 (A) 
^(0) = V6T54Cbcs ^ 20 (A) e ~ 1/150 
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For temperatures 20 K < T < 35 K, the electric field 
E(j) behaves like a power-law only at high current den- 
sities. For smaller values of the current density j, the 
different E(j) curves tend to converge into the T — 35 K 
curve. Further details about this temperature regime will 
be given in the next Section. 

Finally, for temperatures T > 60 K, the E(j) curves 
have a negative curvature in the log E-\osj plot, in agree- 
ment with the interpolation formula ( |l0|) for the bar- 
rier as obtained from the WCPT. The fitting parameter 
U c 1000 K is again consistent with our previous results. 



VI. SUMMARY AND DISCUSSION 

The results obtained in the previous Section for the 
strongly layered Bi2212 single crystal in magnetic fields 
H _L afo-planes are now summarized and further discussed 
within the frame of WCPT. 

For temperatures T < 20 K, the activation barrier for 
vortex motion U(j) depends logarithmically on the cur- 
rent density, while the time relaxation of the current den- 
sity j follows a power-law behavior as given by formula 

ll\ an ap- 



sity 

(|ll|). According to the discussion in Section 
proximately lo gari thmic current dependence in the acti- 
vation barrier ([12]) is obtained within the single vortex 
pinning situation. However, in Section [v| it has been 
shown that for temperatures T < 20 K and small enough 
magnetic fields, the correlation length along the c-axis L c c 
is much smaller than the interlayer distance d. This in- 
dicates that for this regime pinning involves elementary 
pancake- vortices. 

On the other hand, within the most simple approach 
(see Section |n|), for decreasing current densities j the 
activation barrier for single pancakes U (j) is expected to 
be a constant, whereas the measured activation barrier 
is found to be logarithmic up to temperatures T ~ 19 K. 
The non-constant behavior of the measured activation 
barrier U(j) suggests that there are residual interactions 
which were not considered in the most simple approach 
and which lead to an increase of the clastic energy for de- 
creasing j. This argument is also supported by following 
considerations: The collective pinning energy for single 
pancakes,Ej which is the relevant parameter for the deter- 
mination of quantities such as the critical current density 
j c and the depinning energy, is given b y U„ r . ~ e d{j c / j ). 
For the parameters of Table || and |ll| one finds that 



pc 



20 K. Furthermore, the energy which is relevant 
for creep of pancake- vortices is expected to be bigger than 
Up C , but of the same order of magnitude. However, this 
estimated energy is still small as compared to the values 
obtained for the activation energy U(j) plotted in Fig. ||, 
indicating that for creep of pancake-vortices additional 
interactions have to be considered. A possible idea lead- 
ing tp-.coupling of the pancake-vortices into an elastic 
planai3 for decreasing current densities j is the concept 
of variable-range hopping. As discussed in Section H], 



due to the randomness in the energies of the metastable 
state, pancakes will hop over larger distances as the cur- 
rent density j decreases. The shear interaction energy 
CQ^du 2 will therefore grow with increasing hopping dis- 
tance u and for low enough current densities c^du 2 will 
become of the order of the pinning energy U pc . In sum- 
mary, the vortex system is expected to go over from a 
VRH regime (creep of individual pancake-vortices) to a 
2D collective creep regime at low current densities. 
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FIG. 10. Qualitative low field phase diagram of the vor- 
tex state in Bi2212 as proposed in Fig. (The 2D regime is 
not shown in this plot since it is not relevant for the present 
data range). The open circles in the graph are the estimated 
values of the B-field at the center of the sample (B C cntcr) at 
the starting time t = t s ~ 18 • 10 -3 s, just after removing the 
external field H . The vertical arrows qualitatively show the 
time evolution of B CC ntcr during the relaxation of the remanent 
magnetization M rem . The insert represents schematically the 
profiles of the induction B in the remanent state for temper- 
atures between 20 K and 40 K. The shaded areas indicate the 
different pinning regimes throughout the sample (here R is 
the sample radius). In particular, the parameter B13 in the 
insert indicates the value of the induction B in the sample at 
which ID pinning goes over into 3D pinning. The smoothed 
irreversibility line4L in the graph was obtained from the data 
of Schilling et alB 

At higher temperatures (T > 20 K), where bulk pin- 
ning becomes relatively weak, barriers arising from the 
geometry of the sample and from surface effects can play 
an important role. As shown in Section in the low 
magnetic induction regime ("field off" creep measure- 







ments) and for our specimen, these contributions are of 
minor importance as compared to the contributions of 
bulk pinning. 

For temperatures T > 40 K, the activation barrier U (j) 
is found to follow a power-law behavior (Fig. ^|) that can 
be accurately fitted with formula ©. Within WCPT it 
has been shown that a potential of the form of ( |To| ) leads 
to a time relaxation of the current density j as given 
by the interpolation formula (||). From the fit to the 
T = 69 K data in Fig.[|, it is seen that for this temper- 
ature regime the time dependence of the current density 
j is actually well described by formula @ . According to 
WCPT, a power-law potential with an exponent \i ~ 0.6, 
as obtained from the fits, indicates a regime of large bun- 
dle pinning. A list of the obtained fitting parameters is 
given m Table [|. 

No unique functional dependence between U and j 
could be found for the temperature range 20 K < T < 
40 K. Motivated by the present experimental results, the 
very low field regime of the Bi2212 pinning diagram has 
been investigated. Fig. [I] shows the qualitative diagram 
obtained for this regime. For temperatures T above T , 
two different pinning regimes separated by £13 are found 
involving vortex-pancakes and vortex-segments. In or- 
der to describe the relaxation data in this temperature 
regime it is therefore necessary to estimate the induc- 
tion in the center of the sample (-Bcenter) at the time 
t = t s ~ 18T0 _3 s, immediately after removing the exter- 
nal field H. These values are given for all temperatures 
by the empty circles in Fig.|l0|, where the arrows indi- 
cate the time evolution of £? C cntcr during the relaxation 
of the remanent magnetization M rcm . For temperatures 
T < 20 K and T > 40 K, it follows that the whole sample 
is characterized by 0D and by 3D pinning, respectively. 
On the other hand, as shown in the insert of Fig.|l^, for 
temperatures 20 K < T < 40 K, the values of the ini- 
tial field profile in the sample lead to the coexistence of 
two different pinning regimes: 3D pinning in the cen- 
tral region of the sample and ID pinning close to the 
borders. The simultaneous occurrence of two different 
pinning regimes may then provide an explanation why 
no simple functional dependence between U and j could 
be found for this temperature interval. 

From the analysis of the relaxation data we find that 
vortex bundles are strongly pinned against thermal ac- 
tivation. Nevertheless, due to their large size they can 
only sustain low flux density gradients, which means low 
critical current densities j c . On the other hand, vortex 
strings are weakly pinned, but being small in size they 
can sustain relatively high flux density gradients. Much 
higher creep rates are therefore expected for collectively 
pinned vortex lines than for large vortex bundles. Thus, 
a possible interpretation of the vortex dynamics observed 
for the temperature interval 20 K < T < 40 K is the fol- 
lowing: The high relaxation rates which are measured 
at times t < 1 s (see Fig. ||) are mainly the result of the 
strong decay of the flux in the ID regime at the border 
areas of the sample. At times t > 1 s, most of the flux 



has left the sample and only a low flux density gradient 
of vortex strings remains. As discussed in Section III, at 



low current densities j the activated motion of vortex- 
lines involves hops of larger vortex segments by longer 
distances. The weak relaxation rates of the current den- 
sity j measured at times t > 1 s are then explained by 
the growth of the elastic energy for the activated motion 
of vortex strings at low current densities. 

A feature that has recently attracted a lot of interest 
in the literature is the. observation of a second peak in the 
magnetization loop.cll In our Ri2212-sample, as well as in 
several other works on Bi2212,clHl3~Eil this second peak is 
seen for magnetic inductions B of the order of Q /A 2 and 
for tem per atures between 20 K and 40 K. As discussed in 
Section [TV], for magnetic inductions B < <I>o/A 2 the shear 
modulus Cqq starts to decrease exponentially so that, at 
low fields and for temperatures T > 20 K, a ID pinning 
regime can arise (see Fig.[j]). For the ascending branch 
of a magnetization loop it follows that, for fields larger 
than the field of first flux penetration H p , the sample is 
expected to first enter into the ID regime before grad- 
ually going over into 3D. As previously discussed, at a 
fixed temperature T, the value of j c is bigger for the ID 
regime than for the 3D regime. However, since the flux 
in the ID regime has a much faster creep rate as com- 
pared to the 3D regime, it is very important to consider 
the time scale for the measurement of the magnetiza- 
tion loop. For instance, if the time scale were very short 
(t — > 0) , the effects of creep would be negligible and for 
an increasing (decreasing) magnetic field H one would 
expect to measure a decrease (increase) in the magne- 
tization as soon as the magnetic induction B is of the 
order of <E> /A 2 , where the flux in the sample goes over 
from ID to 3D. On the other hand, on the typical time 
scale of the measurement of a magnetization loop, the 
flux in the ID regime is already strongly relaxed while 
the flux in the 3D regime is still close to its configuration 
in the critical state. The value of the magnetization AI 
may then turn out to be smaller in the ID regime than in 
the 3D regime, leading to the characteristic double peak 
in the magnetization loop as measured for temperatures 
between 20 K and 40 K. 

This interpretation is in agreement with previous re- 
ports (Refs. 5l]-|55|) which relate the second peak of the 
magnetization loop to the slower magnetization decay for 
the field range where the peak is observed. Moreover, it 
is in agreement with the results of Ref . [ll] regarding local 
induction measurements on a Bi2212 sample. In the de- 
scending branch of the magnetization loop, at T — 24 K 
and for field values between 330 Oe and 260 Oe, ^uchangc 
of slope of the field profile dB z (x) jdx is observedcil occur- 
ring at various locations inside the crystal starting from 
the edge regions and moving towards the center as the 
applied field is decreased (with B z being the induction 
parallel to the crystallographic c-axis). Within the pre- 
sented low field phase diagram of Bi2212, this change of 
slope is expected to occur at the crossover field B^. 

In conclusion, for temperatures 20 K < T < 40 K, the 
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observation of the second peak in the magnetization loop 
and of the high relaxation rates of the current density j 
for times t < 1 s can both be related to the coexistence 
of two different pinning regimes inside the sample and to 
the strong difference in their relaxation rates. 
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